Thermodynamic Cross-Effects from Dynamical Systems 
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We give a thermodynamically consistent description of simultaneous heat and particle transport, 
as well as of the associated cross-effects, in the framework of a chaotic dynamical system, a general- 
ized multibaker map. Besides the density, a second field with appropriate source terms is included 
in order to mimic, after coarse graining, a spatial temperature distribution and its time evolution. 
A new expression is derived for the irreversible entropy production in a steady state, as the average 
of the growth rate of the relative density, a unique combination of the two fields. 
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The relation between transport processes and chaotic 
models with only a few degrees of freedom (cf. jl]] for re- 
cent reviews) became a subject of active research since 
the rapid progress in dynamical-system theory started in 
the early eighties. First, it was shown that such mod- 
els can faithfully describe particle transport, and be- 
come compatible on the macroscopic level with appro- 
priate macroscopic transport equations Later it was 
found that the irreversible behavior of these processes, 
expressed, e.g., by their entropy production |3|-|6| or fluc- 
tuation relations |?],^| , can properly be obtained in a more 
restricted class of models only — in particular, if one 
wants to keep them low dimensional. Multibaker maps 
§,|0[|-|,|l^||], the extensions of baker maps [0] to a 
macroscopically long array of mutually connected unit 
cells turned out to be an effective tool to understand the 
origin of irreversibility on the level of dynamical systems. 
Besides the possibility of explicit calculations, they lead 
to general findings [Qj^J^l valid also outside the realm of 
multibakers. 

When describing transport in the presence of an exter- 
nal field and/or a density gradient, consistency with the 
thermodynamic entropy balance could only be obtained 
for multibaker maps with a time-reversible, local dissi- 
pation mechanism |^^| (a brief discussion of this notion 
will be given below) . This requirement was interpreted as 
mimicking a thermostatting algorithm (a Gaussian ther- 
mostat, cf. [fl5f ), which is widely applied in nonequilib- 
rium molecular dynamics (NEMD). The entropy balance 
was found to hold for a coarse-grained entropy based on a 
density averaged over regions of small spatial extension. 
A recent paper by Tasaki and Gaspard jl2| shows that 
analogous results can be obtained for area-preserving 
multibaker maps with an energy-dependent phase-space 
volume. This energy, however, is strictly connected to 
the potential of an external field, and not considered as 
an independent driving force. 

In the present paper, our aim is to study transport 
generated by two independent driving forces: density and 
temperature gradients. In addition, we allow for a con- 
stant external field. We are intending to describe a quasi 



one-dimensional system of finite length attached at the 
two ends to different reservoirs, and possibly in thermal 
contact with a thermostat along its extension (cf. Fig. |I|). 
In this general setting we show that thermoelectric phe- 
nomena, i.e., cross effects due to the simultaneous pres- 
ence of two independent driving forces, and the entropy 
balance can properly be modelled by an elementary dy- 
namical system. 
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FIG. 1. Graphical illustration of the transport process con- 
sidered. The system is attached to reservoirs inducing particle 
and heat currents as indicated by the arrow, and along its ex- 
tent heat can be exchanged with a thermostat. 

The multibaker map describes transport along the x 
direction. It represents a deterministic dynamics [the 
(x,p) dynamics] in the single-particle phase space of a 
weakly interacting many-particle system. The cell size a 
partitions the x axis into regions which are sufficiently 
large to characterize the state inside such a cell by ther- 
modynamic state variables and small enough to neglect 
variation of these variables on the length scale of the cells 
(local- equilibrium approximation). Thus, a plays the role 
of a minimum allowed macroscopic resolution. The state 
of the many-particle system is represented by the (par- 
ticle) density g and a so called "kinetic-energy" density 
gT, whose average over cells is related to a local temper- 
ature. For multibaker maps the kinetic-energy density 
is considered as an independent field, i.e., our discussion 
does not rely on the apparance of a momentum conju- 
gated to the x variable. The (x,p) dynamics drives the 
time evolution of the fields, leading to a mesoscopic de- 
scription of the transport process. A possible dependence 
of this dynamics on local thermodynamic averages, and 
the presence of a source term of kinetic energy introduces 
a coupling of the fields. 
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The detailed definition of the model is as follows: The 
multibaker map acts on a domain in the (x,p) plane con- 
sisting of N identical cells labelled by to (Fig. [|). Here, 
x is a position variable, and p is a momentum-like vari- 
able needed to set up a reversible deterministic dynam- 
ics. Every cell has a width a and height 6=1. Af- 
ter every time unit t, every cell is divided into three 
columns (Fig. ||) with respective widths al m , as m and 
ar m fulfilling l m + s m + r m = 1. The right (left) column 
of width ar m (al m ) is uniformly squeezed and stretched 
into a strip of width a and of height l m +i (r m _i) in the 
right (left) neighbouring cell. The middle one preserves 
its area. The map is time reversible in the sense that the 
Jacobian l m +i/r m for a motion from cell to to to + 1 is 
reciprocal to that of the motion from cell m+ltom. The 
(x,p) dynamics is volume preserving when an initial con- 
dition is recovered, but it is contracting on average (since 
motion in the direction of the external field is connected 
with contraction), and hence it is dissipative. Except for 
the coupling to reservoirs at the ends the map is one-to- 
one on its domain. It drives the density g(x,p) and the 
kinetic-energy density g(x,p)T(x,p). Both are advected, 
but — in order to be able to capture a local heating of 
the system — the density gT is also multiplied by a fac- 
tor (1 + rq m ) with q m depending only on the averages in 
cell to and in its neighbors. The dynamics of both densi- 
ties is governed by the Frobenius-Perron equation [jjj of 
the (x,p) map, but after each iteration the gT values in 
cell to are multiplied by (1 + rq m ). The fields g and T 
evolve into fractal distributions whose asymptotic forms 
are described by different invariant measures. In general, 
the width l m , s m , r m of the columns may depend on the 
average values of the fields in the cell and its neighbors 
so that the widths may vary in time and space. 
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FIG. 2. The action of the multibaker map on the coordi- 
nates (x,p) over a time unit r. The values of g(x,p) and 
T(x,p) on the cells and strips [cf. Eqs. ([!]) and (^)] are given 
on the margins for the initial conditions discussed in the text. 



In the spirit of nonequilibrium thermodynamics, we 
also consider the coarse-grained fields g m and T m inter- 
preted as the density and the local temperature of cell 
to, respectively. g m and g m T m are obtained as averages 
of g(x,p) and of the kinetic-energy density g(x,p)T(x,p) 
over cell to. 

We now discuss the time evolution of g and T. For the 
explicit calculation we start with a density g(x,p) and 
a specific kinetic energy T(x,p) taking constant values 
g m and T m in each cell to. This is convenient from a 



technical point of view, and does not lead to a princi- 
pal restriction of the domain of validity of the model as 
discussed in |p|,[tl|. After one step of iteration, the fields 
will be piecewise constant on the strips defined in Fig. ||. 
Due to continuity, the density takes the respective values 
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(1) 



(The prime always indicates quantities evaluated after 
one time step.) Besides the advection by the (x,p) dy- 
namics, the updated values for T on the strips contain 
the source term: 



T ln,r = Trn-l [1 + Tq m ] , 

T Ls = T m [1 + rq m ] , 



(2) 



TL 



rqr, 



In the x variable, the cell-to-cell dynamics of the model 
is equivalent to a random walk with fixed step length a 
and local transition probabilities r m and l m over time 
unit r. Such random walks are characterized by the drift 
v m and the diffusion coefficient D, which stay finite in 
the macroscopic limit a, r — > 0. To be consistent with 
a diffusion type equation for the density, the transition 
probabilities have to scale [|l(| as 

r m -l m = (r/a) v m , r m + l m = (2r/a 2 ) D. (3) 

In order to account for the effect of temperature gradi- 
ents on the local (x,p) dynamics, we allow in the present 
paper for a location dependence of the drift v m , while the 
diffusion coefficient is kept spatially homogeneous. The 
location dependence of the drift is due to a dependence 
on the cell temperature T m and on its discrete gradient: 
v m = v m (T m , (T m+ i - T m )/a). 

The entropy S m of cell to is related to the cell density 
g m . We use the common information-theoretical form 



S m = -ag m \n[g m / g*{T m )] 



(4) 



It is measured in units of Boltzmann's constant, and 
involves a temperature dependent reference density 
g*{T) = T 7 , in close analo gy t o an ideal gas. Here 7 
is a free parameter. In Ref. p3[ it will be demonstrated 
that only this choice of the reference density can be con- 
sistent with thermodynamics. In the present paper, we 
take the form of (||) for granted and explore in how far the 
model leads to proper macroscopic expressions for quan- 
tities appearing in the local entropy balance, i.e., for (i) 
the irreversible entropy production, (ii) the current, (iii) 
the heat and entropy currents, and (iv) the transport 
coefficients. 

Next, we work out the evolution equations of the 
coarse-grained fields and of the entropy S m . Due to O), 
the change of the cell density in a time unit r becomes 



{gj - Q m )/r = -(j m - j m -i)/a, 



(5) 
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where j m = (a/r) (r m g m - l m+ ig m+1 ) denotes the dis- 
crete current density. 

Similarly, an update of the average temperature in the 
cells can be calculated based on the fact that g m T m is a 
kinetic-energy density. Consequently [cf. Eqs. (|l]j2|)], the 
change g' m T' n — Q m T m per time step is obtained as 



Qm^rn QmT„ 



Qm^n 



AqT) _ AqT) 

Jm Jm-1 



, (6) 



where j 



AeT) 



T m jm — Qm+i {T m +\ — T m ) / a (a 2 /r) l m +\ 



is a discrete heat current, and g m T m q m / (1 + Tq m ) repre- 
sents a source of kinetic energy arising from irreversible 
heating. 

In order to find an evolution equation for the entropy, 
we write S' m — S m in the form of a discrete entropy- 
balance equation (cf. ||) 



(7) 



where 



A e S m — S m ^ "5*777 j (8 a ) 

^iSm = {S m ~ ^ ) — {S m — S* m ^), (8b) 

correspond to the entropy flux and the irreversible en- 
tropy production, respectively. Following Refs. j(| and in 
harmony with (4), we define the Gibbs entropy as 

S* (G) = - J dxdp g(x,p) In [g(x,p) T(x,py~<] . (9) 

In view of (8a), the entropy flux is the temporal change 
of the Gibbs entropy. Moreover, (8b) is a meaningful 
measure for the rate of irreversible entropy production, 
since according to the information theoretic interpreta- 
tion of entropies, it describes the increase per unit time 
of the lack of information due to coarse-graining. Hence 

(G) 

it is positive. The dependence of S m — S m on details 
of the coarse graining drops out when taking the time 
derivative [||. 

In view of these considerations, we find the follow- 
ing expressions for the entropy flux and the rate of ir- 
reversible entropy production [cf. Fig. and Eqs. (fOB)] 
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In the last equation we made use of the fact that due to 
the originally homogeneous field distributions, the initial 
coarse-grained and microscopic densities coincide such 

that *S*777 = 5*777 ■ 

A detailed discussion of the entropy current will be 
given elsewhere [ p"3[ . Concerning the entropy production 
in a steady state we only mention the following inter- 
esting relation: The irreversible change of the specific 
entropy per unit time appears to be the average of the 
growth rate of the relative density g(x,p)/g*(x,p), 

i.e., AiS m /(aTg m ) = (cr g / g *) S s, where the growth rate is 
defined as 



t [g(x,p) g*(x,p) 



(11) 



This rule is a direct extension of the result for uniform 
temperature p], where only the growth rate of the den- 
sity appeared. 

Next, we evaluate the entropy balance (for a general 
time-dependent state) in the macroscopic limit. It is ob- 
tained by refining the partition of the finite total length 
L = Na of the chain. Technically this means a, t — > 
with fixed L and diffusion coefficient D. [By this con- 
dition a macroscopic limit v(T,d x T) is obtained for the 
total drift v m ]. The macroscopic spatial coordinate is 
defined as x = am, and the spatial dependence of any 
field -0 = Q or T can be expressed as 
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(10a) 
(10b) 



(12) 



By dividing Eq. (|t]) by ar, we find in this limit a bal- 
ance equation for the entropy density S m /a — > s in the 
form d t s = $ + with 



r (irr) 
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(13a) 
(13b) 



Here, <j>( thermostat ) stands for the entropy flux into the 
thermostat, and j and j( s > are the current and entropy 
current densities, respectively. The former is obtained 
from the time evolution (|5|) of the density, which be- 
comes an advection-diffusion equation in the thermody- 
namic limit with a current j = gv (T, d x T) — D d x g. The 
entropy current can be expressed as 



3 {s) = - [1 + In (QT- 1 )} 3 - 1QD d x T/T, 
and for the flux into the thermostat we obtain 

^(thermostat) = r) gg _ y - 



(14) 



(15) 



The expression for the currents and the irreversible en- 
tropy production are in harmony with thermodynamics 
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Jig . A comparison shows that A = jgD plays the role of 
the heat conductivity, and — T(l + In [pT~ 7 )]) /e is the 
Peltier coefficient (e stands for the unit charge) . The ap- 
propriate form of v (T, d x T) follows from Onsager's reci- 
procity relation JD3]. Thus, we have expressed all the 
kinetic coefficients with system parameters. 

For a system isolated from the thermostat (a non- 
thermostatted system) one has $( thormostat ) = o, which 
fixes q to be vj / A. This choice of q exactly corresponds to 
classical thermodynamics, where the full entropy flux can 
be written as the negative divergence of the entropy cur- 
rent. In this case, the heat leads to a rise of temperature, 
reflected in the source term qT = vjT/X of the evolution 
equation d t T = qT + d x (\ d x T)/(jg) ~ (j/g) 8 X T for the 
local temperature [i.e., the macroscopic limit of @], The 
case q = corresponds to a Gaussian thermostat com- 
monly applied in NEMD simulations Jig ]. Other choices 
of q describe stationary states stabilized by appropriate 
entropy fluxes to (or from) the thermostat. Temperature 
profiles, which are not steady in thermodynamics, appear 
then as steady states by a proper choice of q, implying 
^(thermostat) _^ ^^jj represents a generalization of the 
classical thermodynamic entropy balance. 

To conclude, we point out that to our knowledge trans- 
port with cross effects and the associated entropy balance 
has never been treated in the framework of dynamical 
systems where both the properties of the underlying dy- 
namics and of the thermodynamic time evolution can ex- 
plicitly be worked out (for early efforts based on Gaussian 
thermostats, cf. ]l7p). It is remarkable that a thermody- 
namic description of transport driven by two indepen- 
dent forces could be derived from a model as simple as 
a piecewise-linear map. The minimum requirements on 
the model in order to be consistent with thermodynamics 
are found to be: (i) a time-reversible dissipation mech- 
anism, (ii) inclusion of a new field T describing a local 
kinetic-energy density, (iii) defining the entropy with a 
T-dependent reference density, (iv) inclusion of a source 
term in the time evolution of the field T, which can model 
the effect of Joule's heat. 

Observations which seem to be valid beyond the frame 
of multibaker models are: (a) Dissipation and ther- 
mostatting play different roles. With time-reversible dis- 
sipation we can also describe systems isolated from a 
thermostat (b) Different choices of the source term in 
the temperature dynamics correspond to changes in the 
coupling of the system to the thermostat, (c) A proper 
expression for the total irreversible entropy production 
should involve the steady-state average of the relative 
density g/g*, instead of phase-space contraction, (d) 
In macroscopic steady states, when the coarse-grained 
densities do not change in time, the transition proba- 
bilities and the source term in the dynamics are fixed 
to time-independent values. In this case, one finds a sta- 
tionary entropy balance based on an inhomogeneous low- 
dimensional mapping acting on a spatially extended (i.e., 
macroscopic) domain. 
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